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Evidence of CP violation in the charm sector has been observed re- 
cently by the LHCb and CDF Collaborations. The issue of whether it 
can be accommodated within the standard model (SM) is examined in 
this talk. We conclude that the CP asymmetry difference Aa^fp between 
D° K + K~ and D° tt+tt- is of order -(0.14 ~ 0.15)%. If the im- 
proved theoretical estimate of Aa^ P in the SM remains to be a few per 
mille and the experimental measurement continues to be large with more 
statistics in the future, it will be clear evidence of physics beyond the SM 
in the charm sector. 
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1 Introduction 



In the era of the luminosity frontier, there are basically two strategies for the search 
of New Physics (NP) in the low-energy flavor physics sector: (i) Measure those ob- 
servables which are predicted to be null or almost null in the standard model (SM). 
Examples are CP violation in the decay D~(B~) — > 7r~7r°, lepton number violation 
in t decays, the CP-odd phase (3 S in B s — > J/ifxj), like-sign dimuon asymmetry in 
semileptonic B decays. Although the measurements of these observables will be dif- 
ficult, it will be paid off once the observable of interest is detected. The advantage 
in this strategy is that we don't have to worry much about the SM background, (ii) 
Take a cue from the current anomalies observed at B factories and CDF in the past 
years. The smoking-gun signatures of NP such as the B-CP puzzles of sin 2/3 and 
AAktt (the difference of CP asymmetries in B° — > K~n + and B~ — > K~tt°) and the 
forward-backward asymmetry in the decay B — > K*fi + {i~ are usually of 3 a effects. 
However, many of them were already diminished by the LHCb. Nevertheless, we did 
have very exciting progress in the first strategy, namely, the first evidence of CP vio- 
lation in the charm sector obtained by the LHCb collaboration [1] and corroborated 
subsequently by CDF [2]. 

A nonzero value for the difference between the time-integrated CP asymmetries 
of the decays D° — > K + K~ and D° — > 7r + 7r~ had been reported by the LHCb Col- 
laboration PQ: 



AA CP = A CP (K + K-) - Acp^ti-) = -(0.82 ±0.21 ±0.11)% (LHCb). (1) 



This first evidence of CP violation in the charm sector was later confirmed by the 
CDF Collaboration with the result [2] 



where the direct CP asymmetry, the indirect CP asymmetry, (t) is 

the average decay time in the sample, and r is the D° lifetime. The combination of 
the LHCb, CDF, BaBar and Belle measurements yields c$ P = -(0.025 ± 0.231)% 
and Acicp = —(0.656 ± 0.154)% [3]. The significance of the deviation from zero 
is 4.3 a for direct CP asymmetry. However, for the measurements of AAqp in (TjQ) 
and (T5]), neither LHCb nor CDF has measured the corresponding Acp{K + K~) and 
Acp{t^ + ^) separately. B 

*In 2011, the CDF Collaboration g] has obtained A CP (tt + tt-) = (0.22 ± 0.24 ± 0.11)% and 
A CP {K+K-) = -(0.24 ± 0.22 ± 0.09)%, and hence AA CP = -(0.46 ± 0.31 ± 0.11)% based on a 
data sample corresponding to the integrated luminosity of 5.9 fb _1 . 



AA CP 



(0.62 ±0.21 ±0.10)% 



(CDF) . 



(2) 




(3) 
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It is commonly argued that direct CP violation in singly Cabibbo-suppressed 
(SCS) D decays is of order 0([\V; b V ub \/\V c * s V us \}a s /n) ~ 1(T 4 which is smaller than 
the LHCb measurement by two orders of magnitude. Before claiming new physics 
beyond the SM in the charm sector, it is crucial to have reliable SM estimate of Aa^p. 
In this talk, we would like to address several questions: (i) Can we have a reliable 
theoretical estimate of strong phases and decay amplitudes? (ii) Can the LHCb 
measurement be accommodated within the SM? and (iii) Given the experimental 
results for Aa^ P , can one predict direct CP violation in other SCS charm decays? 



2 Diagrammatic approach 

In order to explore if the surprisingly large CP asymmetries measured by LHCb and 
CDF Collaborations can be explained by the SM, we need to work on a suitable 
theoretical framework. In B physics, there exist several QCD-inspired approaches 
describing the nonleptonic B decays, such as QCD factorization (QCDF), pQCD 
and soft collinear effective theory. However, this is not the case in the D sector. 
Until today we still don't have a satisfactory theoretical framework describing the 
underlying mechanism for exclusive hadronic D decays based on QCD. This is because 
the mass of the charm quark, being of order 1.5 GeV, is not heavy enough to allow 
for a sensible heavy quark expansion. Although from time to time people have tried 
to apply pQCD or QCDF to hadronic charm decays, it does not make too much sense 
to generalize these approaches to charm decays as the l/m c power corrections are so 
large that the heavy quark expansion in l/m c is beyond control. 

Nevertheless, we do have a powerful tool for charm physics which provides a model- 
independent analysis of the charmed meson decays based on symmetry, namely, the 
diagrammatic approach [21 El H]- In this approach, the topological diagrams are 
classified according to the topologies of weak interactions with all strong interaction 
effects included. Based on flavor SU(3) symmetry, this model-independent analysis 
enables us to extract the topological amplitudes and it is complementary to the 
factorization approach. In this diagrammatic scenario, various topological diagrams 
are depicted in Fig. [TJ 

The topological amplitudes T, C, E, A can be extracted from the Cabibbo-favored 
(CF) D — >■ PP decays (in units of 10" 6 GeV) |9]: 

T = 3.14 ± 0.06, C = (2.61 ± 0.08) e - l ( 152±1 )° ; 

E = (1.53l°;° 8 7 ) e ^ 122±2 )°, A = (0.39^) e* 31 ^" . (4) 

The topological amplitudes C and E extracted from the data are much larger than 
those expected from naive factorization. Color suppression in the C amplitude and 
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(c) P, Pg w (d) S, ftw (a) PE, PE av (h) PA, PA„ 

Figure 1: Topology of possible flavor diagrams: (a) color-allowed tree T, (b) color- 
suppressed tree C, (c) QCD-penguin P, (d) singlet QCD-penguin S diagrams with 2 (3) 
gluon lines for M<i being a pseudoscalar meson P (a vector meson V), (e) ^-exchange E, 
(f) VF-annihilation A, (g) QCD-penguin exchange PE, and (h) QCD-penguin annihilation 
PA diagrams. The color-suppressed EW-penguin -P^w an< ^ color-favored EW-penguin Pew 
diagrams are obtained by replacing the gluon line from (c) and all the gluon lines from (d), 
respectively, by a single Z-boson or photon line. The EW-penguin exchange -PEew and 
EW-penguin annihilation -PAew diagrams are obtained from (g) and (h), respectively, by 
replacing the left gluon line by a single Z-boson or photon line. 



helicity suppression in E are alleviated by nonfactorizable effects [j and final-state 
rescattering. Recall that the ^-exchange amplitude scales as l/m c and vanishes 
in the limit of m c — > oo. Hence, the sizable E extracted from the data is clear 
evidence for the importance of l/m c power corrections. Considering the existence of 
an abundant spectrum of resonances at energies close to the mass of the charmed 
meson, we can even make quantitative and qualitative statements on FSI effects in 
weak annihilation. As shown in Refs. [101 [11], the effect of resonance-induced FSI's 
can be described in a model-independent manner in terms of the mass and width 
of the nearby resonances. It is found that the E and A amplitudes are modified by 
resonant FSI's as (see, e.g., Ref. [TT] ) 

E = e + (e 2 ^ - 1) (e + |) , A = a + (a 2 * - 1) (a + ^) , (5) 

with 



mo — m R + iT/2 

T Nonfactorizable contributions will render the 1/N C term in the factorizable amplitude C sup- 
pressed. 
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where the short-distance contributions to the jy-exchange amplitude E and W- 
annihilation amplitude A are denoted by e and a, respectively. Therefore, even if 
the short-distance weak annihilation is turned off, a long-distance jy-exchange (W- 
annihilation) contribution can still be induced from the tree amplitude T (C) via FSI 
rescattering in resonance formation. 



3 Singly Cabibbo-suppressed decays 

As far as direct CP violation is concerned, it occurs only in SCS D decays. But why? 
If the hadronic decay of charmed mesons proceeds only through tree diagrams, then 
CP violation will not occur as only the quarks of the first two generations will get 
involved. In order to have CP nonconservation in the SM, all three-generation quarks 
should participate in the weak decay process. The heavy b quark enters the D decays 
through the loop diagram. Since the penguin diagram does not occur in the CF and 
doubly Cabibbo-suppressed (DCS) decays, direct CP violation at tree and loop levels 
will manifest only in the SCS decays. Now writing the SCS decay amplitude as 



Amp = A rf (tree + P d ) + A s (tree' + P s 



(7) 



with \ p = V* V up , we obtain the general expression of CP asymmetry as 



n dir ~ 

a CP ~ 



2Im(A d A:) 
IA.I 2 



D 


sin 5 = 2 






D 






sin 7 




f 




f 



sin 5 = 1.2 x 10" 



sin 5, 



where T is the dominant tree amplitude and D can be a tree or penguin amplitude, 
7 and 5 are the weak and strong phases, respectively. For D being a tree amplitude, 



one will expect 10 > a CP > 10 . For the interference between tree and penguin 
amplitudes, it is expected that a^°p P ^ ~ 10~ 4 as \P/T\ ~ 0(a s (m c )/7c) provided that 
the relative strong phase is close to maximal. 



3.1 SU(3) symmetry breaking 

In order to have a reliable calculation of important task is to understand the 

SU(3) breaking effects in SCS D — > PP decays. A most noticeable example of SU(3) 
breaking lies in the decays D° — > K + K~ and D° — > 7r + 7r~. Experimentally, the rate 
of D° — y K + K~ is larger than that of D° — > 7r + 7r~ by a factor of 2.8 [12J, while they 
should be the same in the SU(3) limit. This is a long-standing puzzle since SU(3) 
symmetry is expected to be broken merely at the level of 30%. 

The conventional wisdom for solving the above-mentioned long-standing puzzle is 
that the overall apparently large SU(3) symmetry violation arises from the accumu- 
lation of several small and nominal SU(3) breaking effects in the tree amplitudes T 
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and P (see e.g.,[13\). From the recent measurements of AAcp by LHCb and CDF, 
we learn that penguin diagrams in SCS decay channels do play a crucial role for 
CP violation. This leads to the conjecture that penguins may also explain the rate 
disparity between D° — > K + K~ and D° — > 7r + 7r~. 

To discuss SU(3) breaking effects, we express the decay amplitudes of D° — > 7r + 7r~ 
and K + K~ in terms of AP and XP: 

A (D° 7T+7T-) = X d {T + E + P d + PE d + PA d ) m + X s (P s + PE s + PA s )^ 

= ^(X d -X s )(T + E + AP) 7Tn -h b (T + E + ZP) n7T , (9) 



where 



AP = (P d + PE d + PA d )-(P s + PE s + PA s ), 

SP ee (P d + PE d + PA d ) + (P s + PE s + PA s ) . (10) 



Likewise, 



A(P° A+p-) = i(A. - A d )(T + P - AP)^ - ^A b (T + P + SP) JCA . . (11) 

As far as the rate is concerned, we can neglect the term. Now we consider 
three different cases: (i) SU(3) symmetry holds for the T and P amplitudes. Then a 
sizable AP with a negative real part is needed to contribute constructively to K + K~ 
and destructively to 7r + 7r~; that is, AP ~ ^Te~ i200 ° [H]. If [/-symmetry breaking in 
the amplitudes T + P follows the pattern [15J 



(T + EU = (T + P)(l + ~e T ), (T + P) KK = (T + P)(l — ie T ), (12) 



where e T is a complex parameter with \e T \ G (0,0.3), it has been shown in [15] 
that the relation |AP/T| ~ 0.5 still holds roughly, (ii) The realistic symmetry 
breaking in T and P amplitudes does not necessarily follow the pattern given in 
Eq. ( fl2l . Indeed, according to the factorization approach, SU(3) violation due 
to decay constants, meson masses and form factors leads to the robust relation 
T^/Ttttt = 1.32ai(PP)/ai(7r7r) w 1.32. Assuming SU(3) symmetry again for the 
P amplitudes, we will have |AP/T| ~ 0.15. This is also the case considered in [16J. 
(hi) The large rate difference between D° — > K + K~ and D° — > 7t + tt~ is entirely 
accounted for by SU(3) violation in the T and P amplitudes and AP is negligibly 
small. Owing to the observation of D° — > K°7C° through ^-exchange and penguin 
annihilation diagrams and the smallness of AP theoretically, we have argued in [H] 
that the last scenario is preferred. The W -exchange amplitudes can be fixed from the 
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following four modes: K + K , 7i + n , 7r°7r° and K°K° . Neglecting AP and terms, 
a fit to the data yields two possible solutions 

(I) E d = 1.19 e il5 -°°£, E s = 0.58 e~ iU - 7 ° E , 
(II) E d = 1.19e il5 -°°£, E s = l.Q2e~ i9 - 8 °E , (13) 

where E g refers to the H^-exchange amplitude associated with cu — >■ gg (g = (i, s). 
The corresponding x 2 vanishes as these two solutions can be obtained exactly. 



4 Direct CP violation 
4.1 Tree- level CP violation 

Most theory papers after the LHCb measurement focused on the decays D° — > K + K~ 
and D° — > tt + tt~ and emphasized the importance of interference between tree and 
penguin amplitudes for CP violation. However, direct CP asymmetry can occur at 
the tree level in many of the SCS decays of D mesons in addition to the loop-level 
CP violation [17]. The strong point of the topological approach is that the magnitude 
and the relative strong phase of each individual topological tree amplitude in charm 
decays can be extracted from the data [18]. Hence, it becomes possible to make a 
reliable estimate of a^ r . . Larger tree- level CP asymmetries can be achieved in those 
decay modes with interference between T and C or C and E. For example, aj^ ee ' ) is 
of order (0.7 - 0.8) x 10" 3 for D° 7r°r/ and D+ ->• K + r] (see Tabled]). 
Direct tree-level CP violation in D° — > K°TC° is given by 

for the two solutions for E d and E s given in Eq. ( IT3l) . For comparison, a d ^ e \D° — > 
K°K°) = 1.11 x 10~ 3 is predicted in [19]. 

From Table [T], we see that almost all the predicted tree-level CP asymmetries in 
[19] are of opposite signs to ours. This is ascribed to the phase of the IV-exchange 
amplitude. For CF D — > PP decays, its phase is (122 ±2)° with \ 2 — 0.29 per degree 
of freedom [Eq. (jl}]. For SCS decays, the phases of E d and E s [see Eq. ( TIB"]) ] lie in the 
range of 107° ~ 137°. Therefore, the W -exchange amplitude in this work is always 
in the second quadrant, while the E amplitude in [19] lies in the third quadrant from 
a global fit to all the data of 28 CF and SCS D — > PP branching fractions with 
X 2 = 7.3 per degree of freedom. As a result, the imaginary part of E in [T9] has a 
sign opposite to ours, and this explains the sign difference between this work and 

r (tree) 
for °>dir ■ 
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Decay Mode 


a^ 00 ' (this work) 


a£T\l9 


(tot) 

a dir ( this work ) 




Expt. 


D u - 


> 7T + 7T 








0.95 ±0.04 


0.68 


2.0 ±2.2 


D°- 


>7r°7r° 








0.80 ±0.04 


0.20 


1±48 


D°- 


» ir°r] 


0.82 ±0.03 


-0.33 


0.08 ±0.04 


-0.55 




D°- 


> ttV 


-0.39 ± 0.02 


0.54 


0.01 ± 0.02 


1.99 




D°- 


> 7777 


-0.28 ± 0.01 


0.28 


-0.58 ± 0.02 


0.08 








-0.42 ± 0.02 


0.28 


-0.74 ± 0.02 


0.08 




D°- 


> 7777' 


0.49 ± 0.02 


-0.30 


0.54 ± 0.02 


-0.98 








r\ 1 i"» no 

0.38 ± 0.02 


—0.30 


0.34 ± 0.02 


—0.98 




D°- 


> X+X- 


U 


U 


— U.4z ± U.U1 


— U.oU 


n 1 1 n 
— Z.O ±1.1 






n 
u 


n 
U 


n c^q 4- n no 
— U.oo ± u.Uz 


— U.OU 




D°- 


>K°K° 


-0.73 


1.11 


-0.63 ± 0.01 


1.37 








-1.73 


1.11 


-1.81 ± 0.01 


1.37 




D+ - 
















29 ± 29 


D+ - 


-> 7T + ?i 


0.35 ±0.06 


-0.54 


-0.74 ±0.06 


-0.52 


17.4 ± 11.5 


D+ - 


-> 7T+t/ 


-0.21 ±0.04 


0.39 


0.33 ±0.07 


1.52 


-1.2 ± 11.3 


D+ - 




-0.07 ±0.06 


-0.14 


-0.39 ±0.04 


-1.00 


-1.0 ±5.9 


D+- 




0.07 ±0.06 


0.14 


0.45 ±0.03 


1.00 


66 ±24 


Df- 




0.01 ±0.11 


0.33 


0.94 ±0.10 


0.72 


266 ± 228 


Df- 


-> #+?7 


-0.71 ±0.05 


-0.19 


-0.61 ±0.05 


0.83 


93 ± 152 


D+- 




0.36 ±0.04 


-0.41 


-0.28 ±0.12 


-1.78 


60 ± 189 



Table 1: Direct CP asymmetries (in units of 10 3 ) of D — > PP decays. The first (second) 
entry in D° — > rjrj, rjrj' , K + K~ and K°K° is for Solution I (II) of E d and E s [Eq. fjl3H ] . For 
QCD-penguin exchange PE, we assume that it is similar to the topological E amplitude 
[see Eq. (|18j) ]. World averages of experimental measurements are taken from Ref. [3J. For 
comparison, the predicted results of ci^* and a^*° ™ P2] are a ^ so presented. 



4.2 Penguin- induced CP violation 



Direct CP violation does not occur at the tree level in D° — > K + K~ and D° — > 7t + tt~ 
decays. CP asymmetry in these two modes arises from the interference between tree 
and penguin amplitudes denoted by aj^r ■ Specifically, we have 



(t+p) 

h dir ' 



7T 7T 



1.2 x 1(T 3 

-1.2 x icr 3 



P s + PE 8 + PA S 



T+E+AP 
P d + PE d + PA 



sin 5- 



7T7T 1 



T+E-AP 



sin 6 



KK 



KK ' 



(15) 



where 5 nn is the strong phase of (P s + PE S + PAg)^ relative to (T + E + AP)^ 
and likewise for Therefore, in SU(3) limit we have the relation (K + K~) = 

The QCD penguin contributions can be estimated in QCD-inspired approaches 
such as QCD factorization (QCDF) and pQCD or in the factorization approach. In 
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Decay Mode 


(tree) 


u dir 


(t+ 


pa) 


_(tot) 
a dir 


D° -> 7T + p~ 





0. 


.08 


-0 


.60 


-0.52 


D° -» 7T-p+ 





— u. 


Oct 


-0 


.22 


-0.28 


D° -> 7T p 





— U. 


no 


-0 


.74 


-0.76 







— 0. 


no 
.08 





.60 


0.52 







u. 


.0 < 





.22 


0.29 


no v fr-oiz* 


73 


0. 


.73 


n 
\j 


73 


73 


nO __> t>-*0 


— U. lo 


— 0. 


.73 


— U 


. 1 3 


— U. / 3 


D° -> n°co 





-0. 


.01 





.53 


0.52 


D° -> ttV 


















D° ->■ 


0.19 


0. 


.19 





.50 


0.50 


D° r/u; 


-1.07 


-1. 


.05 


-0 


.91 


-0.89 




















L>° Tip 


-0.53 


-0. 


.55 


-0 


.22 


-0.24 


d° ip° 


0.59 


0, 


.59 





.21 


0.21 



Table 2: Same as Table [T] except for D — > PV decays where the superscript (t + p) denotes 
tree plus QCD penguin amplitudes, (t + pa) for tree plus weak penguin annihilation (PE 
and PA) amplitudes and "tot" for the total amplitude. Due to the lack of information 
on the topological amplitudes Ap and Ay, no prediction is attempted for D + — > PV and 
D+ ->■ PV decays. 



P '' 0.24e-* 154 °, 0.010, ' 



T \ T 

± / 7T7T \ -L 

P d ^> --■ il52° (Pd — Ps\ n nnn „-»35 



QCDF we find 



0.24e-' iM , % = 0.009e"^. (16) 

^ ' KK \ 1 / KK 

The difference in the d- and s-loop penguin contractions turns out to be very small 
compared to the tree amplitude. It follows that 

.35e* 170 °, (W) =0.24e 4l70 °. (17) 



T + EJm \t + e 

Hence, 5^ « 5 KK = 170°. This leads to a { ^ p \n+n-) = 6.7xl0- 5 and a^ p) (K + K~) = 
— 4.9xl0 -5 . Therefore, QCD-penguin induced CP asymmetries in D° — > 7r + 7r _ , K + K~ 
are small mainly due to the almost trivial strong phases 5 nn and 8kk- 

For power corrections to QCD penguins, we shall consider QCD-penguin exchange 
PE and QCD-penguin annihilation PA. At the short- distance level, weak penguin 
annihilation contributions are found to be smaller than QCD penguin with the hier- 
archy \P\ > \PE\ > \PA\. For example, (PE/T)^ ~ 0.04 and (PA/T)^ 0.02. 
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1J 




u K~ u 



Figure 2: Long-distance contribution to D° — > n + n~ through a resonant-like final- 
state rescattering. It has the same topology as the QCD-penguin exchange topological 
diagram PE. 

As for long-distance contributions to weak penguin annihilation, it was pointed out 
in [18] that a SCS decay, for example, D° — > 7r + 7r~, could proceed through the weak 
decay D° — > K + K~ followed by a resonant-like final-state rescattering as depicted in 
Fig. [2j It has the same topology as the QCD-penguin exchange topological graph PE. 
Just as the weak annihilation topologies E and A, it is expected that weak penguin 
annihilation will receive sizable long-distance contributions from FSI's as well. Hence, 
we shall assume that PE, PEp and PEy are of the same order of magnitude as E, 
E P and Ey, respectively; namely, 

PE « E, PE P « E P , PE V « E v . (18) 

The calculated direct CP asymmetries for various SCS D — > PP and D — > PV 
decays are summarized in Tables Q] and [2} respectively. We conclude that the direct 
CP asymmetry difference between D° — > K + K~ and D° — > tt + tt~ is about — (0.139 ± 
0.004)% and —(0.151 ± 0.004)% for the two solutions of IV-exchange amplitudes, 
respectively. A similar prediction of —0.118% was also obtained in [19]. Since in the 
SM, Aa,Qp arises mainly from weak penguin annihilation, we can vary the amplitude 
PE to see how much enhancement we can gain. Even with the maximal magnitude 
\PE\ ~ T and a maximal strong phase relative to T, we will get Aa*p = -0.27% 
which is still more than 2a away from the current world average. 

4.3 Extraction of penguin amplitudes from the diagrammatic 
approach? 

The cautious reader may notice that while we have relied on the topological tree ampli- 
tudes extracted from the data to perform the analysis, we apply the theoretical model 
to compute the penguin amplitudes. The question is then can we extract penguin 
amplitudes from SCS decays within the framework of the diagrammatic approach? 
As shown in [16], in principle the penguin amplitudes J2 P =d,s^p(Pp + PEp + PAp) 
can be extracted from D° — > 7i + 7i~,7i iT ,K + K~ decays, while J2 P =d,s \{Pp + PEp) 
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can be inferred from D + — > K + K and Df — > tt + K , ir K + provided that the tree 
amplitudes are those given in Eq. (TjJ. Hence, J2 P =d,s \PAp can be determined. How- 
ever, unlike topological tree amplitudes, it will be difficult to extract the topological 
penguin amplitudes reliably from the data. This is because the use of the topological 
approach relies heavily on SU(3) symmetry which leads to negligible penguin ampli- 
tudes in D decays. Consequently, the extraction of penguin amplitudes depends on 
SU(3) breaking effects in tree amplitudes. Indeed, we have shown in Sec. 3.1 that the 
difference in penguin contractions AP is sensitive to how SU(3) symmetry breaking 
is treated in the tree amplitudes T and E . Therefore, we shall rely on theory for the 
estimation of the penguin contribution. 

5 Attempts for the SM interpretation 

Even long before the LHCb experiment, it had been advocated in [20] that the ob- 
servation of CP violation in D decays would not necessarily be a signal of NP. If 
certain hadronic matrix elements are enhanced, in analogy to the AI = 1/2 rule of 
kaon decays, then CP violation will be observable in strangeness conserving decays. 
This has motivated many authors to conjecture that the large direct CP asymmetry 
difference Aa^p observed by LHCb and CDF could be explained or marginally ac- 
commodated in the SM [13 [TBI Ell E21 1231 EH- Most of them argued that the penguin 
matrix elements were enhanced by some nonperturbative effects or unforeseen QCD 
effects. The idea is that the QCD penguin amplitude is of A J = 1/2 transition and it 
will be enhanced relative to the AI = 3/2 amplitude through the A J = 1/2 rule en- 
hancement. However, we shall show below that the enhancement of the A$ amplitude 
relative to A 2 in charm decays arise predominately from the tree amplitudes. 
Combining the amplitudes 



with the decay amplitude of D° — > n + n given in Eq. ([9]), we obtain the isospin 
amplitudes 




A 







)=[\ d {2T -C + 3E)+Y1 3A P (P P + PE P + PA P )}, 




(20) 
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Experimentally, the ratio I-A0/A2I is determined from the relation 

(21) 



t(d+ 


-> 7T+7T ) 


3 


l^2| 2 


T(D° -> 7T + 7T- 


) + r(£>° -)• 7r°7r°) 


~ 2 


|A)| 2 + |^ 2 | 2 



It follows that |^4o/^2| is equal to 2.5 ± 0.1 in charm decays and 22.4 ± 0.1 in the 
kaon system. To a good approximation, 



An 



Ai 



J_ 2T - g + 3P + 3AP 

V2 T + C • (22j 



Neglecting the penguin contribution for the moment and assuming SU(3) symmetry 
for the topological amplitudes so that we can apply Eq. (jJJ) to obtain |A /A 2 | = 3.8 . 
When SU(3) symmetry breaking in T and C amplitudes |T4] are taken into account, 
one would have 



A 



A, 



1 1.92T - 0.79(7 + 3E + 3AP 
71 0.96T + 0.79C ^ ^ 



and the ratio will be reduced to 3.2 . At any rate, it is evident that the AI = 1/2 
enhancement over the AI = 3/2 one in D decays arises predominately from the tree 
amplitudes. As the predicted ratio tends to be too large compared to experiment, this 
means that the difference AP in penguin contractions should contribute destructively 
to |A /A 2 |. This is very different from the kaon case where the predicted ratio due 
to tree amplitudes is too small compared to experiment, and it is quite obvious that 
a large enhancement of the penguin matrix element is needed in order to explain the 
discrepancy between theory and experiment for the ratio |^4o/^4.2| in K — > ttti decays. 

Under the assumption of large enhancement of the penguin amplitude |SP [see 
Eq. ( 1T0|) ] relative to the tree one by a factor of 1/e', Brod et al. argued that SP could 
explain Aa^p, while the difference AP explained the large disparity in the rates of 
D° — > K + K~ and D° — > 7t + tc~ [15]. This would require that AP be of the same 
order as the tree amplitude. A consistent picture for large penguins in SCS D — > PP 
decays is emerged provided that the relevant amplitudes scale as 

T~0(1), ^SP~0(l/e'), AP ~ 0(e u /e') ~ 0(1) , (24) 

with ejj being a [/-spin breaking parameter. While this scenario sounds appealing, we 
notice that the realistic symmetry breaking effects in amplitudes T and E does not 
follow the pattern depicted in Eq. ( fT2l) . For example, the factorization approach leads 
to the robust result {T^/T^l w 1.32 . Assuming SU(3) symmetry for E amplitudes, 
we will have |AP/T| ~ 0(0.15) rather than |AP/T| ~ 0(1). And it will become 
smaller when symmetry breaking in E is taken into account. In realistic model 
calculations, |AP/T| is very small and negligible. The quantity AP comes from the 
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difference of d- and s-loop contractions induced mainly from the 4-quark operator 0\. 
Penguin contraction amplitudes read 



I A ( 



CfCVs 

AnN r 



Cl 



4 777, 2 

3 ln 7T + 3 ~ GM 2( r P' fc2 ) 



+ ■ ■ 



where p = d, s, r» = m 2 /m 2 cl k 2 is the momentum squared carried by the virtual gluon 
and the perturbative loop function G(r, k 2 ) is given by 



G(r, k 2 



-4 / duuil — u) In 
Jo 



r — u(l — u)- 



mi 



(25) 



Hence, 



pd 
^4,6 



ps 
^4,6 



1.6 



1,6 



[GK/m2,fc 2 )-GK/m^P)]. 



(26) 



By varying /c 2 between m 2 and m 2 /4, one can see that AP is rather small, as shown 
in Eq. ([TBI). 

It is often stated in the literature that due to the large l/m c corrections, it is 
reasonable to assume an enhancement of the hadronic matrix elements by a factor 
of, say, 3 ~ 5. However, we notice that this enhancement is not applicable to the 
color-allowed tree amplitude T as its prediction based on the factorization approach 
agrees with experiment even before enhancement. For the amplitude E, its large l/m c 
power corrections arise from final-state rescattering. Hence, it is the decay amplitude 
rather than the hadronic matrix element that gets enhanced by l/m c effects because 
the long-distance FSI cannot be expressed as a single hadronic matrix element of local 
4-quark operators. 

We digress here to make a side remark on the so-called AA^-k puzzle related to 
the difference of CP asymmetries of B~ — > K~tt° and B° — > K~n + . The decay 
amplitudes of B — > Kit in terms of topological diagrams read 



A(B~ ->■ K-tt°) 



P + T + -P£ W + P A , 

_L( P + T + C + P EW + ^ew + A + P A ). (27) 



We notice that if C, Pew and A are negligible compared with T, it is clear from 
Eq. (j27|) that the decay amplitudes of K~it° and K~tt + will be the same, apart from 
a trivial factor of l/Vz. Hence, one will expect that Acp(K~tt°) « Acp(K~n + ), 
while they differ by 5.6a experimentally, AA KlT = Acp(K~ti ) — A C p(K^n + ) = 
0.124±0.022 [3J. A large penguin P cannot explain the AA^^ puzzle as it contributes 



equally to both B~ ->■ K'tc and B° -> K~ 



71 



This puzzle will be resolved provided 
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that c/T is of order 1.3 ~ 1.4 with a large negative phase, where c = C + Pem 
[\c/T\ ~ 0.9 in the standard short- distance effective Hamiltonian approach) [25J. 
There exist two popular scenarios for achieving a large complex c: either a large 
complex C or a large complex Pew (see [25J for references and details). However, 
for a large complex Pew one needs NP beyond the SM because Pew is essentially 
real in the SM as it does not carry a nontrivial strong phase. In principle, one 
cannot discriminate between these two scenarios in penguin-dominated decays as 
C and Pew are comparable in magnitude due to the CKM enhancement for the 
latter. Nevertheless, these two scenarios will lead to very distinct predictions for tree- 
dominated decays where Pew ^ C. For example, the decay rates of B° — y tt°tt°, p°7r° 
will be substantially enhanced for a large C but remain intact for a large Pew- Owing 
to the large branching fractions observed for B° — y 7r°7r°, p°7r°, it is more appealing 
to resolve the B —y Kir CP puzzle using a large complex C instead of invoking NP 
on electroweak penguins. 



6 New Physics Effects 

If the experimental measurement continues to be large with more statistics in the 
future, it will imply new physics beyond the SM in the charm sector. Then it will be 
important to explore possible NP scenarios responsible for such large direct CP asym- 
metries. 

The surprising LHCb measurement has inspired many different analyses based on 
a variety of NP models. Scenarios for NP effects at tree level include flavor-changing 
coupling of a SM Z boson ETJ, flavor-changing neutral currents induced by a 
leptophobic massive Z' boson [2B1 [27], two Higgs-doublet model [27], color-singlet 
scalar model [29] , color-sextet scalar model [27J [30] , color-octet scalar model [27] and 
fourth generation model [311 [23]. Models with NP in QCD penguins at the loop 
level have been constructed as well, including new fermion and scalar fields [27] and 
the chirally enhanced chromomagnetic dipole operator [SB]. CP violation induced by 
supersymmetric P-parity violating interactions has also been considered in |32j . 

It is known that NP models are highly constrained by D°-D mixing, 
mixing and CP violation in the kaon system characterized by the parameter e'/e [33J. 
Many of the tree-level NP models are either ruled out or in tension with flavor-related 
experiments [27]. As pointed out in [26J, a large NP contribution to the AC = 1 
chromomagnetic dipole operator is probably the best candidate to explain the LHCb 
and CDF results as it is least constrained by all current data in flavor physics. 

In [13] we have considered two possibilities of new physics effects, namely, large 
penguins and large chromomagnetic dipole operator, and studied their phenomeno- 
logical consequences in the CP asymmetries of SCS D —y hh decays, seeing if there 
are discernible differences in the two scenarios. 
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6.1 Large penguins 

It is known that a large penguin of order 3T can explain the observed Aa^p [T5| [16]. 
More precisely, the penguin amplitudes 

1 / 2.9Te i85 ° for Solution I 

2 \ 3.2Te i85 ° for Solution II 1 6) 

with nearly maximal strong phase can accommodate the measurement of direct 
CP asymmetry difference between D° — > K + K~ and 7r + 7r~ (see Table [3]). This 
can be easily seen from the approximated relation 



Aa d j P « -2.4 x 10" 



T + E 



sin 5, (29) 



derived from Eq. (II 5p . where we have applied the approximation of a^p (7r + 7r~) rs 
—a < ^ I P (K + K~). Using the above large penguin |EP as input, the predicted direct 
CP asymmetries for other SCS charm decays are summarized in the second column of 
Table |3j We see that many modes, such as — > tt + K°, tc°K + , K + r]', are expected 
to yield direct CP asymmetries of a similar magnitude, at a few per mille level. 

Penguin amplitudes could be enhanced due to some nonperturbative effects in 
the SM 1 or NP effects beyond the SM. The predictions in Table |3] are presented 
irrespective of the origin of enhancement. However, whether or not the large penguins 
are subject to the constraints from D°-D mixing, K°-K mixing, etc. remains to be 
investigated. 



6.2 Large chromomagnetic dipole operator 

The authors in [26J have argued that a large chromomagnetic dipole operator could 
be the best NP candidate to explain the data while satisfying most flavor physics con- 
straints at the same time. Although the chromomagnetic dipole operator Og g is sup- 
pressed by the charm Yukawa coupling, the hadronic matrix element (MxM2\0& g \D) 
scales as m^/h 2 , where k 2 is the square of momentum transfer of the gluon and is 
of order m 2 . As a consequence, the matrix element is independent of m c ; that is, it 
is enhanced by a factor of l/m c . After absorbing the factor m c into the definition 
of 0$ g from new physics loop contribution, such as a low-energy supersymmetry sce- 
nario discussed in detail in [26J , the corresponding Wilson coefficient looks like being 
enhanced by a factor of v/m c , where v is the vacuum expectation value of the Higgs 

* Authors of [16] have introduced an additional phenomenological penguin amplitude Pb in order 
to accommodate the measured Aa^p. However, there is no need to make this assumption as the 
penguin amplitude can be recast to XdPd + A S P S = ^(A^ — X s ){Pd — P s ) — \Xb{Pd + P s )- The second 
term on the RHS of the above relation is the so-called Pb in [TB] . 
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Decay Mode Large penguins Large c.d.o. 
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(0) 


D+ 


— > 7l + 1] 


-3 


24 


(-3.62) 


-5 


35 


(-3.67) 


D+ 


— > n + t]' 


2 


97 


(3.34) 





93 


(0.59) 


D+ 




-2 


95 


(-3.28) 





37 


(0.29) 


Df 


->• TT+fT 


3 


29 


(3.66) 


-0 


47 


(-0.35) 


Dt 




4 


57 


(5.08) 


4 


40 


(3.14) 


Dt 


K+T] 


-0 


58 


(-0.57) 


1 


59 


(0.94) 


Dt 


—> K + T]' 


-5 


16 


(-5.79) 


1 


76 


(1.39) 



Table 3: Direct CP asymmetries (in units of 10 3 ) of SCS D — > PP decays estimated in 
the scenarios with large penguin contributions and large chromomagnetic dipole operator. 
The parameters SP and c^ p are chosen to fit the data of Aa^ p : ±£P = 2.9Te iS5 ° and 
c np = o.017e il4 ° for Solution I, ±£P = 3.2Te i85 ° and c^ p = 0.012e il4 ° for Solution II. 



field and used to represent the typical new physics scale. On the contrary, the D°-D 
mixing induced by 0$ g is suppressed by a factor of m 2 c /v 2 . This illustrates why the 
dipole operator can escape the constraint from D°-TJ° mixing. In short, we need NP 
to enhance the Wilson coefficient cs g and to induce a sizable imaginary part. This can 
be realized in the supersymmetric models where the gluino-squark loop contributes 
a major part of c$ g [21], the disoriented A terms and split families are the sources of 
flavor violation [26J, or the flavor structure of the trilinear scalar couplings is related 
to the structure of the Yukawa couplings via approximate flavor symmetries [34] , or 
the supersymmetric realization of partial compositeness [SB] . 

For the purpose of illustration, we shall take c^ 5 P ~ 0.012e* 14 ° which fits to the 
data of Aa,Q P . The calculated CP asymmetries for the other modes are listed in the 
last column of Table |3j It is interesting to notice that while a large chromomagnetic 
dipole operator leads to large direct CP asymmetry for D° — > 7r°7r°, 7c r], the predicted 
CP violation for D° -> n°r]', D+ -> n+rf, K+K° and D+ -»■ n+K°, K+rj' is much 
smaller than that in the large penguin scenario. Therefore, measurements of the 
CP asymmetries of the above-mentioned modes will enable us to discriminate between 
the two different NP scenarios. 
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6.3 Discrimination between SM and NP interpretations 

In view of two possible explanations of large CP asymmetry difference Aa^ P by the 
SM and by NP, it is natural to ask if there are experimental tests that allow us to 
distinguish between the SM and NP interpretations. There are a few proposals for 
this purpose: (i) Suppose NP has AI = 3/2 contributions. If electroweak penguin 
and isospin breaking effects are negligible, then a measurement of nonzero a^p in 
D + — > 7T + 7r° will be a signal for AI = 3/2 new physics [36]. Some sets of isospin sum 
rules can be constructed along this direction, (ii) If a large chromomagnetic dipole 
operator is responsible for Aa^ P , the corresponding electromagnetic dipole operator 
will also get enhanced. Direct CP asymmetries in radiative D — > P + P~7 decays with 
the invariant mass Mpp close to the p or the peak can be achieved at the level of 
several percent [37]. Hence, it was claimed in [37] that evidence of |app 7 | > 3% would 
be a clear signal of physics beyond the SM and a clear indication of new CP -violating 
dynamics associated to dipole operators. 

Since the SM interpretation of large Aa^ p relies on the large penguin amplitude 
which is about 3 times bigger than the color-allowed tree amplitude, the predictions 
shown in Table |3] will also serve the purpose of discriminating the SM from the NP 
model based on the chromomagnetic dipole operator. As elaborated on in the previous 
subsection, SM and NP in dipole operators can be distinguished as their predictions 
for CP asymmetries in D° — > 7r°7r , 7t°t], • ■ ■ , — > K + i]' appear to be "orthogonal" 
to each other. 

6.4 Other decay modes of interest 

Since the large direct CP asymmetry difference in D° — > K + K~ and D° — > tt + tt~ 
has been observed, one may ask if such an effect also occurs in other decay modes. 
From Table |3] and the measured branching fractions of SCS D — > PP decays [12], it 
turns out that D + — > K + K° and — > n + K° will be good targets for this purpose. 
However, direct CP violation in these charged D decays will be contaminated by 
K°-K mixing due to the neutral kaon K$ in the final state. The decays D° — > 
D + — > 7r + ?/ ' and Df — > K +r q ( - ) will be the next of interest. Observation of 
CP violation in the decay D° — > K°K° will be quite interesting, but its rate is just 
too small. 

As for SCS D — > VP decays, one should measure CP violation in D° — > n + p~ , n~p + , 
K + K*~ , K~K* + to see if their CP asymmetry differences are comparable to that 
measured in D° — > K + K~ and D° — > 7r + 7r~. It is worth mentioning that AAqp is 
measured to be (0.51 ± 0.28 ± 0.05)% in D + (pn + by Belle [38]. Naively, its direct 
CP violation is expected to vanish as D + — > 07r + proceeds only through the color- 
suppressed tree diagram Cp where the subscript P indicates that the pseudoscalar 
meson contains the spectator quark of the charmed meson. Nevertheless, this de- 
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cay also receives singlet QCD-penguin contribution S [see Fig. [Tfd)] governed by the 
effective parameters 03 and 05. The interference between Cp and S will allow the 
presence of CP asymmetry in this mode. 

From the experimental point of view, it is important to search for CP -violating 
effects in three- and four-body charm decays. For example, the Dalitz-plot analysis 
of D + — > K + K~tt + allows us to differentiate strong phases from the weak ones across 
the Dalitz plot. However, theoretically it will be difficult to estimate these effects 
due to the lack of information on l/m c power corrections. The model-independent 
diagrammatic approach has not been proved to be applicable to three-body decays. 

7 Conclusions 

Based on the topological diagram approach for tree amplitudes and QCD factorization 
for a crude estimation of perturbative penguin amplitudes, we have studied direct 
CP asymmetries in singly Cabibbo-suppressed D decays within the framework of the 
SM and concluded that the CP asymmetry difference Aa^p between D° — > K + K~ 
and D° — > 7r + 7r~ is of order —(0.14 ~ 0.15)%. We have explored the phenomenological 
implications of two new physics scenarios for explaining the observed CP asymmetry 
in the charm sector, one with large penguin amplitudes and the other with a large 
chromomagnetic dipole operator. 
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